Utilization of photon orbital angular momentum in the low-frequency radio domain 
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We show numerically that vector antenna arrays can generate radio beams that exhibit spin and orbital an- 
gular momentum characteristics similar to those of helical Laguerre-Gauss laser beams in paraxial optics. For 
low frequencies (< 1 GHz), digital techniques can be used to coherently measure the instantaneous, local field 
vectors and to manipulate them in software. This enables new types of experiments that go beyond what is pos- 
sible in optics. It allows information-rich radio astronomy and paves the way for novel wireless communication 
concepts. 
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Classical electrodynamics exhibits a rich set of symmetries 
[1], and to each Lie symmetry there corresponds a conserved 
quantity [2, 3]. Commonly utilized conserved electromag- 
netic (em) quantities are the energy and linear momentum, 
where the underlying symmetries, under Poincare transforma- 
tions, are homogeneity in time and in space, respectively. 

Another conserved quantity, manifesting the isotropy of 
space, is the em angular momentum, whose mechanical prop- 
erties were predicted theoretically in 1909 [4] and demon- 
strated experimentally in 1936 [5]. A collection of nonrela- 
tivistic, spinless, classical particles with linear momenta pj^^^^ 
has angular momentum J^^^^^ _ £.(x; — xq) x pf^^^^, where 
xo is the moment point. When this system interacts with 
em radiation with spin angular momentum (SAM) S^"^ and 
em orbital angular momentum (0AM) L^"^, the total angu- 
lar momentum, J^^^ = J^^^^^ ^ jem^ ^^^^^ jem ^ gem ^ L^"^ = 

£o / d^x (x — Xo) X (E X B), is conserved [6]. Hence, for a 
fixed S^"^, a change in J^^^^^ ^jjj result in an opposite change 
in L^"^, observed as a rotational (azimuthal) Doppler shift 
[7] . This shift is distinct from the translational Doppler shift, 
which is a manifestation of the conservation of linear momen- 
tum. As shown in Ref. 8, using an angular momentum flux 
representation [9], it is always possible to separate a beam 
into S^"^, which depends only on the local polarization struc- 
ture, and L^"^, which depends on the gradient of the fields; see 
also Ref. 10. 

During the past few decades the use of em 0AM (beam vor- 
ticity) has come to the fore in optics [11] and in atomic and 
molecular physics [12]. However, while SAM (wave polariza- 
tion), generated by proper phasing of the two legs in a crossed 
dipole or by using helix antennas, has been used routinely for 
at least half a century, 0AM has not yet been utilized to any 
significant degree in radio physics [13] or its applications such 
as radio astronomy [14]. The use of radio 0AM is currently 
being contemplated for detection of ultrahigh energy neutri- 



nos interacting with the Moon [15], studies of radio wave in- 
teractions with the atmosphere and ionosphere [16-18], and 
radar probing of the Sun [19, 20]. 

We propose to use antenna arrays for generating and de- 
tecting both SAM and 0AM in radio beams and show numer- 
ically how this works. In such arrays one needs access to the 
complete 3D vectors of the radio em field over an area which 
is large enough to intersect a substantial fraction of the radio 
beam. This requires vector antennas, e.g., tripoles [21, 22]; 
crossed dipole antennas will be useful for beam directions 
nearly perpendicular to the antenna planes. Using digital sam- 
plers connected directly to each vector antenna, the local, in- 
stantaneous 3D radio field vectors themselves can be mea- 
sured coherently up to the GHz range, enabling their manipu- 
lation, including J^"^ processing, entirely in software. This 
is in contrast to infrared and optical frequencies for which 
current detectors are incapable of measuring first-order field 
quantities. There phase-coherent down-conversion to the low- 
frequency radio domain might provide a solution. 

Lasers often use Laguerre-Gaussian (LG) modes [23] in 
which the phases of the electric and magnetic vector fields in 
a plane perpendicular to the beam axis have an l(p dependence 
where / is an integer and (p is the azimuthal angle. This means 
that for / 7^ the phase fronts of LG modes are not planar but 
helical. As shown in Ref. 24, this implies that LG beams carry 
an 0AM of Ih per photon. In the paraxial approximation, the 
LG modes form a complete basis set for light beams [23]. 

A pure 0AM state radio beam of frequency (0 and energy H 
has a beam axis 0AM component Lf^ = IH/co and the fields 
have an azimuthal phase dependence of exp(i/(p), where / is 
an integer as in an LG beam. In order to study the possibility 
of using 0AM in radio, we consider an antenna array and as- 
sume for simplicity that each antenna is located equidistantly 
along the perimeters of circles (see Fig. 1). The antennas are 
fed the same signal, but successively delayed relative to each 
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FIG. 1: (color online). Radiation patterns for radio beams generated 
by one circle of 8 antennas and radius A plus a concentric circle with 
16 antennas and radius 2A; all antennas are 0.25 A over the ground. 
Notice the influence of / on the radiation pattern. Here / = (upper 
left), / = 1 (upper right), / = 2 (lower left), and / = 4 (lower right). 



FIG. 2: Samples of instantaneous electric field vectors E across the 
main lobes of the beams in Fig. 1 (same / values and plotting order). 
The size of an arrow is linearly proportional to the local |E|. As 
expected for 0AM carrying beams, the phase of the em field changes 
by 1271 for a full turn around the beam axis. 



other such that after a full turn around the antenna array axis, 
the phase has been incremented by 1271. The far field inten- 
sity patterns, calculated with the software package NEC2 [25] 
which solves Maxwell's equations for a given set of antenna 
currents, are displayed in Figs. 1,3, and 4 and are very similar 
to those obtained in paraxial optics. Fig. 2 displays the instan- 
taneous E field vectors across the main lobes for different /. 
Theory predicts that j = l^s= (OJf^/H [24] . Table I lists the 
results obtained for / = 0, 1 , 2, 3 and 5- = — 1 for a radio beam 
generated by an antenna array. The agreement is excellent. 

The superposition of two coaxial 0AM states was modeled 
numerically by using two concentric antenna rings with dif- 
ferent radii. In the inner ring we increase the overall phase 
between successive antennas by 27tl\/n\ . In the outer ring we 
use the phase increment 27r/2/^2- The variables n\ and n2 are 
the number of antennas in each ring, respectively, and l\ and 
I2 are the phase increment factors corresponding to the 0AM 
number / in the LG beams. The resulting antenna radiation 
patterns are shown in Fig. 3 for three different combinations 
of l\ and h. These patterns would be difficult to synthesize 
without resort to the 0AM technique. For comparison, inten- 
sity patterns for 0AM carrying LG beams for the same l\ and 
I2 combinations are also shown. 

In the general (linear) case, the total 0AM in a beam is a su- 
perposition of several 0AM states. This superposition can be 
decomposed into pure 0AM states via a discrete Fourier trans- 
form. In particular, one has to integrate the complex field vec- 
tor weighted with exp(— i/(p) along a circle around the beam 
axis. Since there will be only a finite number of antennas 
along the integration path, there is an upper limit on the largest 



0AM number that can be resolved. Namely, |/| <#{antennas 
on a circle around the beam axis}/2. We have assumed that 
the beam axis is centered on the antenna array, but this is not 
necessarily true. The beam axis is determined by the emitter 
only, not by the receiving antenna array. Therefore an incom- 
ing radio beam might not overlap the antenna array perfectly 
so that only an asymmetric spatial part of the beam can be 
analyzed. A disadvantage in this case is that we do not have 
exact information of the phase along the whole circle around 
the beam axis. We have to extrapolate from accurate mea- 
surements of the field vectors within the finite size antenna 
array to the field vectors around the beam axis. This extrap- 
olation yields an uncertainty A/ > RA(p/D^ {R:^ D), where 
A 9 is the smallest phase difference in radians that can be re- 
solved, D is the diameter of the antenna array, and R is the 
distance from the beam axis to the antenna array. Since / 
is integer valued, the uncertainty does not matter as long as 
it is less than 1/2. While losing information about the indi- 



TABLE I: Scaling of coJ^'^/H as a function of / for a right-hand 
circular polarized beam (5 = — 1) formed by a ring array of 10 crossed 
dipoles. Array radius D = X, antennas 0.1 A over perfect ground, 
polar angle = 0. 
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FIG. 3: (color online). Beams obtained by superimposing two dif- 
ferent 0AM states. The upper three panels show the radiation pat- 
terns for the antenna array, and the lower three panels show the 
corresponding intensity patterns, head on, calculated for Laguerre- 
Gaussian beams. The leftmost are for h = \ and I2 = 2, the middle 
ones are for h = I and I2 = 4, and the rightmost are for li =2 and 
/2 = 4. Notice the good agreement between the patterns obtained 
with the antenna array model and the paraxial LG beam model. 

vidual 0AM states, it is possible to measure estimated larger 
0AM numbers up to |/| <#{ antennas along a circle segment 
of length DjnR/D, (7? > D). 

The Poynting vector of a radio beam with 0AM has a he- 
lical phase structure and spirals around the main beam axis 
with a pitch angle a/ = arctan(A//2;r7?). This angle can be 
resolved if / can be resolved, if A/ < 1/2. In this case one 
observes multiple images of a single point source where each 
image corresponds to a pure 0AM state /. If / cannot be re- 
solved, the multiple images blend together yielding a smeared 
spot. Conversely, if the individual 0AM states are resolved, 
the smeared spots are resolved as multiple images (one for 
each /), enabling self-calibration techniques that sharpen the 
radio image via the use of 0AM. 

Inserting typical values of D = 100 km, Acp = 2n x 
17360°, into Al > RA(p/D, and requiring that Al < 1 /2, one 
obtains Rmax = 3000 km for which the 0AM can be resolved. 
It should therefore be possible to probe 0AM processes in 
the Earth's ionosphere and lower magnetosphere and also to 
actively induce plasma vorticity and strong toroidal electric 
currents there [16]. In particular, a radio beam that interacts 
with a turbulent medium will carry information on the vortic- 
ity of this medium allowing for remote radio imaging of the 
turbulence [17]. But, it is not likely that radio beams from 
distant astronomical objects are so narrowly focused that they 
can be probed for 0AM unless the array extends into space. 
On the other hand, if the ionosphere adds 0AM to a beam 
from a distant source, this added 0AM can be measured and 
compensated for. 

Should, however, the radio emitting object rotate fast and 
have sharp discrete lines in its emission spectrum, the angular 
momentum of the emitted em waves can be measured indi- 
rectly via shifts and splittings of the spectral lines. The shifts 
result from the rotational Doppler effect (0 — COq = {I -\- s)Q.^\ 




FIG. 4: (color online). Radiation patterns for a circularly polarized 
beam propagating obliquely (0 = 25°) with / = (left) and / = 1 
(right), generated by phasing the individual elements of a ten-tripole 
array in free space. This illustrates that with a tripole array it is pos- 
sible to control electronically both the beam direction and /. This is 
not possible with arrays of single or crossed dipoles. Note that for 
/ ^ 0, there will be an on-beam-axis minimum which can be useful 
to block out a bright object when observing faint surrounding objects 
[26], e.g., in the solar corona [20]. 



where ^ly is the projection of the rotation frequency onto the 
wave vector k [7]. Decomposing into pure spin states, the 
discrete emission spectrum will be decomposed into one for 
^ = +1 and one for ^ = — 1. These two spectra should almost 
coincide in their spectral lines, except for an overall shift. The 
relative overall shift between the two spectra is equal to twice 
the rotational frequency of the emitter. Once has been read 
off, one can search for spectral lines that are separated exactly 
by (and integer multiples thereof). Each of these spectral 
lines corresponds to a specific 0AM state. 

When pointing the beam in a direction other than the z axis 
(orthogonal to the plane of the dipoles), vector sensing an- 
tennas such as tripoles are much preferred. They have the 
advantage that the three orthogonal antenna currents can be 
viewed as the x, y, and z components of a composite antenna 
current vector j(^, x) which, under software control, can be ro- 
tated into any given direction in space, in the receive case even 
after the fact. To generate off-zenith em beams with non-zero 
/ values, the tripole is indispensable since one cannot control 
the direction of the beam with linear dipoles or crossed dipoles 
and still maintain the desired Icp phase dependence; see Fig. 4. 
This is because the phase difference between the elements is 
used to point the beam. The only way to overcome this limita- 
tion would be to mechanically rotate the antenna array. Such 
an arrangement would be infeasible. 

In summary, we have shown by theoretical and numeri- 
cal modeling that a number of the photon SAM and 0AM 
signatures observed in paraxial optics can be generated in 
the radio frequency range where, for low enough frequencies 
(< 1 GHz), modern digital radio techniques can be used to 
measure coherently the 3D electric (and magnetic) vectors in 
the beams. This opens up new kinds of fundamental electro- 
magnetic radiation experiments. Furthermore, minor changes 
of the design of the large low-frequency multi-array radio tele- 
scopes coming on line (LOFAR [27]) or in the planning stage 
(SKA [28]), would enable them to utiHze not only SAM but 
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also 0AM, which would increase their resolution, sensitiv- 
ity, interference tolerance, and overall usefulness. The LOIS 
(LOFAR Outrigger in Scandinavia) Test Station in southern 
Sweden has already implemented the radio 0AM techniques 
and is currently being set up for the first proof-of-concept ex- 
periments [19]. 

Since information can be encoded as 0AM states [29] that 
span a much larger state space than the two- state SAM space 
[30], radio 0AM techniques hold promise for the development 
of novel information-rich radar and wireless communication 
concepts and methodologies. Furthermore, it is conceivable 
that signatures related to radio OAM or other conserved quan- 
tities in Dirac's symmetrized form of the Maxwell-Lorentz 
em theory might provide clues on the existence of magnetic 
monopoles [31]. Finally, it should be mentioned that invari- 
ants are very useful for assessing the stability and robustness 
of numerical simulation codes. 
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